Population change map

library(testmap)

Use the interactive map to check population changes.

points_bau4326 <- points_bau |> 
  #filter(BASIN_N == "Orinoco") |>
  arrange(BASIN_N, SUBBASI, BB_ID, REACH_ID) |> 
  filter(row_number() %% 5 == 1) |> st_as_sf() |> st_transform(4326)
points_bau4326_low <- points_bau4326
# colour palette
pal <- colorFactor(
  palette = c("#7274C1", "#A3720E"),
  domain = points_bau4326$flag_50_42y
)
 
# interactive map
leaflet::leaflet(data = points_bau4326) |> 
  addTiles() |>
  addCircleMarkers(color = ~pal(flag_50_42y), 
             stroke = FALSE, fillOpacity = 0.4, 
clusterOptions = markerClusterOptions(), 
group = "points_bau4326"
) |>                        
  groupOptions("points_bau4326", zoomLevels = 1:8) |> 
  addCircles(data = points_bau4326_low, 
             color = ~pal(flag_50_42y), 
group = "points_bau4326_low"
) |>                        
  groupOptions("points_bau4326_low", zoomLevels = 9:15)